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Charmed Meson Decays and QCD Sum Rules 
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The current status of the QCD sum rule predictions for charmed mesons is overviewed. 



1 Introduction 

Weak decays of charmed mesons provide a useful testing 
ground for nonperturbative methods, such as QCD sum 
rules [ 1 1. Contrary to B decays, in D decays the hadronic 
matrix elements are multiplied by experimentally known 
CKM parameters: V cs is extracted from charm-tagged W- 
decays and V c d from neutrino-nucleon production of charm 
[|2). Therefore, the QCD sum rule predictions for the 
hadronic parameters of D decays can be directly compared 
with experimental data. This comparison allows one to 
gain more confidence in the sum rule results for the B- 
decay hadronic matrix elements that are used to extract the 
poorly known CKM parameters. 

The outline of the QCD sum rule method and the predic- 
tions for B decays are discussed in the context of CKM 
physics in [ 3 1 (in sect. 4.2 of chapter III and in sects. 2.1.3, 
2.1.4 of chapter IV), see also [@j. Here I will concentrate 
on the corresponding results for the charm sector. A de- 
tailed review of QCD sum rules can be found in [0. 

2 Determination of the c-quark mass 

One of the first applications of QCD sum rules was the es- 
timate of the charmed quark mass. According to the orig- 
inal method l"|6|. one employs the correlator of two cy^c 
currents. Due to dispersion relation, the n-th derivative of 
the correlator in q 2 (the momentum transfer squared) is re- 
lated to the n-th power moment of the hadronic e + e~ — ► cc 
cross-section a c : 




The moments M„(q 2 ) are calculated in QCD in a certain 
range of {n,q 2 }, far off-shell (q 2 < 0), in terms of the vir- 
tual c-quark loops, taking into account the c-quark inter- 
actions with perturbative gluons and nonperturbative gluon 
condensate [ 1 1. The results depend on m c , a s and the gluon 
condensate density. The hadronic cross section in Eq. Q 
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is saturated by the charmonium resonances J/ y/, , and 
at y/s > 2mo by the open charm-anticharm production. 

In Table 1 the recent results for the MS mass m c (m c ) ex- 
tracted from the sum rule moments Q with the 0(a 2 ) ac- 
curacy are presented and compared with the predictions of 
other methods. 



Table 1. 



m c (m c ) (GeV) 


Ref. 


method 


1.304(27) 
1.275(15) 


[Q 
[1 


QCD SR 


1.230(90) 


[ED 


QCD SR+NRQCD 


1.190(110) 


[Q3 




1.370(90) 


[ED 


FESR 


1.210(70)(65)(45) 


[El 


rriB — mo + HQET 


1.300(40) (200) 


[El 


lattice QCD (average) 



Despite a reasonable agreement between all four QCD SR 
predictions within theoretical uncertainties, one has to keep 
in mind that the quoted estimates of m c (m c ) are obtained 
using M„ with different {n,q 2 } and making different as- 
sumptions. In [ 7 1, following the original analysis of [[6] 
the lowest moment n = 1 at q 2 = is selected, having a 
little sensitivity to nonperturbative effects but demanding 
an accurate knowledge of the cross section G c (s) above the 
open charm threshold. In the analysis of [|8), n ~ 10 and 
q 2 < are used, so that the gluon condensate contribu- 
tions become important. Finally, in rl9l 1101 . for the mo- 
ments with large n, an attempt is made to account for the 
resummed Coulomb effects that are not accessible in the 
relativistic calculation of Yl(q 2 ). An ansatz for the spec- 
tral density of the correlator is adopted, combining the full 
QCD answer with the resummed NRQCD spectral density 
at large and small c-quark velocities, respectively. This 
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choice emphasizes Coulomb versus gluon-condensate ef- 
fects in a relatively light cc quarkonium, an issue which de- 
serves further studies (for a critical discussion see [0). In 
order to further improve the sum rule determinations of m c , 
one needs more precise measurements of the J /\j/,y/' 
leptonic widths and of the open charm cross section in 
e + e~ . 

Furthermore, let me remind that the vector charmonium 
channel provides a possibility to check the quark-hadron 
duality approximation, the key element of QCD sum rules. 
Replacing in the integrand in Eq. Q the hadronic cross 
section by the calculated spectral density of the correlator 
at s > so, one still successfully fits the moments, even if so 
is shifted close to the open charm threshold, indicating that 
the "semi-local" quark-hadron duality also works in this 
channel. 

3 f D from SVZ Sum Rules 

One of the great advantages of the sum rule method is a 
possibility to easily switch from c to b quark and vice versa, 
in the analytical expressions. An important example is the 
calculation of the Z)-meson decay constant fo defined as 
f D m 2 D = m c (0\diy 5 c\D + ). The SVZ sum rule is derived 
from the correlator of two ciy^d currents, applying duality 
approximation for excited D states and Borel transforma- 
tion. The result has the following schematic form: 

n ""' x I ;(y \i \ 

E E V ) Snk{m c ,M 2 , s °,n c )L n (n c ) (2) 

n=0,3,... \*=0,1,.. % ) 

where S„k are the calculable short-distance coefficients. In 
particular So* are given by perturbative heavy-light loops 
in order (Lq =1); L„>3 are the universal long-distance 
parameters (vacuum condensate densities) with dimension 
n. In the above, M ~ m c A (where A = mo — m c ), is the 
Borel parameter characterizing the average virtuality of the 
c-quark in the correlator, Sq is the quark-hadron duality 
threshold , and fl c ~ M is the factorization scale. The fact 
that L„>3 ~ (Agco)" an d Snk ~ (1 /M)" allows one to retain 
a finite numbers of terms in the sum over n ( n max = 6 is 
already providing a sufficient accuracy). 

The fo — > /b transition in the sum rule (|2) is realized by 
replacing 

m c -> m b , m D -> m B , -> s%, fa -> \l h (3) 

with the scale-dependence given by the relevant 
renormalization-group factors. The c <-> b correspon- 
dence does not mean however that both sum rules for fg 
and fu are equally accurate. The numerical hierarchies of 
corrections are different in the appropriate ranges of Borel 



parameters, so that generally the D-meson sum rules are 
less stable. Nevertheless, since fs is of a crucial impor- 
tance for Z?-physics it is very useful to check the method by 
comparing the sum rule prediction for fu with experiment. 

In Table 2 two recent sum rule determinations of fo are 
presented (for a more detailed overview including older re- 
sults see [ 5 1). Comparison with the lattice QCD results re- 
veals an encouraging agreement. In [Q4| the perturbative 
part is taken into account in 0((X S ), and the quoted theoret- 
ical uncertainty is largely determined by the c-quark mass 
interval. The one-loop pole mass rn c , "" p = 1.3 ±0.1 GeV 
was taken, which overlaps with the lower part of the m ] c !oop 
interval obtained from the values of m c (m c ) given in Ta- 
ble 1 . One possibility to improve the sum rule is to use the 
recently obtained 0(a£ ) results for the heavy-light correla- 
tor [[15 1, providing the coefficient S02 in Eq. (|2). This was 
done in [ 16 1 in the framework of HQET In full QCD so 
far only fs was recalculated [ 17 1 with the O(aj) accuracy. 

In addition, to have more confidence in the power expan- 
sion of the correlator, it would be useful to calculate the 
d = 7 correction proportional to a combination of quark 
and gluon condensates. A better determination of m c , in- 
clusion of 0(a s ) 2 corrections in full QCD, together with 
a systematical use of m c (m c ), are the remaining resources 
of improvement for the sum rule result for fo- More diffi- 
cult is to assess the "systematic" uncertainty related to the 
quark-hadron duality approximation in the D meson chan- 
nel. 



Table 2. 



fo (MeV) 


Ref. 


Method 


200(20) 


[El 


SVZ, 0{a s ) 


195(20) 


[\m 


SVZ+HQET, 0(a 2 ) 


203(14) 




Lattice QCD average (quench.) 


226(15) 


[M 


(unquench.) 



Having obtained a prediction for fo one is not yet able to 
compare it with an experimental number, because only fo s 
is measured, the latest result [ 19 1 being : 

f Ds = 285 ± 19 ± 40 MeV. (4) 

Importantly, QCD sum rules also predict the fo, I fo ratio 
in terms of m s and (ss) /(qq), (q = u,d) the ratio of strange 
and nonstrange quark condensates. The rather old results 
collected in [ 5 1 yield an interval: 

f D JfD = 1.11 -1.27, (5) 

to be compared with the recent averages fojfo = 
1.12(2) [1.12(4)] of lattice quenched [unquenched] QCD 
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r II 81 with a smaller uncertainty (see also [ 20 1). Multiply- 
ing the sum rule prediction from Table 2: fo = 200 ±20 
MeV by the ratio 101 we obtain: f Ds = 240 ±40 MeV, in 
the ballpark of the experimental interval @. Improving the 
latter and measuring fo will provide more decisive checks. 



4 D — > n,K form factors from LCSR 



Measuring the semileptonic D — > nl\>i decay distribution 



with a DK = 0.01 



+0.11 

-.07 ' 



2N3/2 



dq 2 



247T 3 



\.f£M 2 )] 2 

+0(mj), 



(6) 



at m? < q 2 < (mo — m K ) 2 and dividing out \V ce j\ one is 
able to reproduce the experimental values of the D — > n 
form factor f^ n (q 2 ). With a sufficiently large statistics, 
of semileptonic decays with / = ji, the scalar form factor 
fo n (q 2 ) entering the chirally suppressed 0(mj) part of the 
decay distribution can also be extracted. The Cabibbo en- 
hanced D^f Klvi decays provide D—^K form factors with 
an even better accuracy. 

Heavy-to-light form factors are calculated using QCD 
light-cone sum rules (LCSR). One of the first applications 
of this sum rule technique was the calculation of DK 
in [|^. The updated results for fo % oK including higher 
twist terms [ 22 1 and 0(a s ) corrections H23I can be found 
in N14I . The form factor f^ >n was calculated in H24I . 

The sum rules for B — > n are obtained by the same replace- 
ment (fjjl. Let me emphasize that this transition is done 
from one finite mass to the other. Contrary to HQET rela- 
tions between B and D form factors, no heavy-quark mass 
approximation is involved. Since /g_ (q 2 ) is used to extract 
\V„b\ from B — > Klv, a comparison of the sum rule predic- 
tions for D — ► % form factors with experimental data will 
ensure more confidence in the LCSR method. 

The most recent LCSR result f DK (0) = 0.65 ±0.11 ob- 
tained in [ 14 1 takes into account the twist-2 term in NLO 
and twist-3,4 contributions in LO (in the expansion of the 
underlying vacuum-pion correlator in the pion distribution 
amplitudes with growing twist). From the same correla- 
tor, using double dispersion relation one has access to the 
D*Dn coupling go*Dn [ 1221 1251 . predicting the product 
fD*fDgD*DK- Using the SVZ sum rule result for fo quoted 
above, the D* -pole contribution to the D^n form factor is 
calculated. The two predictions of LCSR are used t il 41 to 
fit the D — > K form factor in the whole kinematical region 



< q < (mo — m n 
dispersion relation: 



to the simple ansatz [ 26 1 inspired by 



fi 



065 ±0.11 



D% 



(l-q 2 /m 2 D , ) ( 1 - a DK q 2 /m 2 D , ) 



(7) 



Interestingly, the sum rule results 
are consistent with the D*-pole dominance for the form 
factor. The predicted integrated decay width F(D° — > 
n-l+Vi)/\V cd \ 2 = 0.13 ±0.05 ps" 1 is in a reasonable 
agreement with the experimental number [ [2) r(D° — > 
7t-e + V e )/\V cd \ 2 =0.174 ±0.032 ps" 1 . At the zero mo- 
mentum transfer the LCSR prediction also agrees with the 
recent lattice result [|27) / D k(0) = 0.57(6)±{] ^. On the 
other hand, the lattice calculation [ 1271 suggests that the 
contribution of excited D* states is not small. To assess the 
reliability of these theoretical predictions, one has to wait 
until the D — > nl\>i decay distribution is accurately mea- 
sured. 

Finally, the LCSR result fmD^K form factor [ [14| 
is /dk(0) = 0.78 ±0.11 at m. ( (lGeV) = 150 MeV, very 
sensitive to the strange quark mass. The correspond- 
ing integrated width T(D° -> K~ l+Vi)/\V C!i \ 2 = 0.094± 
0.036 ps _1 , is in a good agreement with experiment [|2) 
r(D°^/ir/+V/)/|V c . v | 2 =0.087 ±0.004 ps -'. 

Concluding this section, let me briefly discuss the impor- 
tant issue of the D*Dn coupling (see also [ 28 1). The re- 
cent first measurement of the total width of D* by CLEO 
collaboration [ 29 1: T tut (D*) = 96±4±22 keV yields 
for this coupling go*Dn — 17.9 ± 0.3 ± 1.9 (using the def- 
inition of Ref. [ |22|). The LCSR prediction \ l22l l25l. 
gD-DK = 10 ±3.5, is obtained by dividing the calculated 
product fofD*gD*DK by the two decay constants, fo and 
fo* , extracted from 2-point SVZ sum rules. Taking into ac- 
count the estimated theoretical uncertainty, the upper limit 
of the LCSR prediction is go*D% = 13.5, still 25% lower 
than the central value of the CLEO number. Meanwhile, 
the first lattice QCD prediction go*Dn = 18.8±2.3+2o [ 
30 1 agrees with the CLEO result. If the future measure- 
ments (although extremely difficult !) and lattice calcula- 
tions confirm the large value of this coupling, one has to 
clarify the status of the LCSR prediction. Having in mind 
that all other sum rules discussed above use one-variable 
dispersion relations, one might suspect that certain com- 
plications arise in the double dispersion relation used in 
LCSR for the D*Dn coupling. More specifically, the sim- 
plest quark-hadron duality ansatz (one resonance plus con- 
tinuum) in both D and D* channels may be too crude. One 
possible scenario was recently discussed in [ 3 1 1 : assuming 
a partial cancellation between the contributions of excited 
and ground D,D* states in the dispersion relation, one gets 
an increase of the LCSR coupling. Without going into fur- 
ther details, let me only mention that the magnitude of the 
coupling go<Dn and the shape of the form factor f^ n (q 2 ) 
are closely related. Suppose the form factor is dominated 
by the D*-pole contribution: 



Dm 



2m D *(l -q 2 /m Dt ) 



(8) 



Taking for go*Dn the CLEO central value and multiply- 
ing it with fo* = 250 MeV (within the lattice QCD pre- 
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diction [ 32 1) one obtains a semileptonic width F(D° — ■> 
n~l + Vi)/\V c d\ 2 = 0.37 ps _1 which is about two times 
larger than the experimental width quoted above [ |2) . To 
make the strong coupling measured by CLEO consistent 
with the semileptonic width, one needs a substantial nega- 
tive interference between the D*-pole and excited D* states 
in the dispersion relation for the form factor (as also no- 
ticed in [ 1311 ). resulting in a visible deviation of the D — > 
TilVi decay distribution from the D*-pole dominance. 

5 Rare D decays 

Exclusive rare D decays such as D — > ,2y, etc. will 

become important highlights in the future high-statistics 
charm physics experiments. Being suppressed in the Stan- 
dard Model, these decays are promising indicators of new 
physics. The long-distance amplitudes of rare D decays 
in the Standard Model still lack a QCD based analysis. I 
believe, QCD sum rules in both two-pont (SVZ) and light- 
cone versions can essentially help in solving this problem, 
One example is the LCSR prediction for weak radiative de- 
cays obtained in f 1331 : 

BR{D+ -> p + y) = 2.7 -10- 6 , 
BR(D° ^ p°j) = 3.0 -10- 6 , 

BR(D S -> p+y) = 2.8 -lO" 5 , (9) 

where all numbers have an (9(50%) accuracy. This analysis 
can be further improved and extended to the other rare D 
decay channels. 
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